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ABSTRACT 

We present new Suzaku and Fermi data, and re-analyzed archival hard X-ray data from 
INTEGRAL and Sw//f-BAT survey, to investigate the physical properties of the luminous, 
high-redshift, hard X-ray selected blazar IGR J225 17+22 17, through the modelization of its 
broad band spectral energy distribution (SED) in two different activity states. Through the 
analysis of the new Suzaku data and the flux selected data from archival hard X-ray obser- 
vations, we build the source SED in two different states, one for the newly discovered flare 
occurred in 2005 and one for the following quiescent period. Both SEDs are strongly domi- 
nated by the high energy hump peaked at 1 20 - 1 22 Hz, that is at least two orders of magnitude 
higher than the low energy (synchrotron) one at 10 n -10 14 Hz, and varies by a factor of 10 
between the two states. In both states the high energy hump is modeled as inverse Compton 
emission between relativistic electrons and seed photons produced externally to the jet, while 
the synchrotron self-Compton component is found to be negligible. In our model the observed 
variability can be accounted for by a variation of the total number of emitting electrons, and 
by a dissipation region radius changing from within to outside the broad line region as the lu- 
minosity increases. In its flaring activity, IGR J22517+2217 shows one of the most powerful 
jet among the population of extreme, hard X-ray selected, high redshift blazar observed so 
far. 

Key words: Galaxies: active - quasars: individual (IGR J22517+2217) - radiation mecha- 
nisms: general - X-ray: galaxies 



1 INTRODUCTION 

In the past years, Integral-IBIS (Ubertini et al. 2003) and Swift- 
BAT (Barthelmy et al. 2005) catalogs of hard X-ray selected extra- 
galactic sources have been published (Sambruna et al. 2007; Bird 
et al. 2010), opening a new window in the study of blazars. 

The hard-X-ray selection produces a set of diverse sources, in- 
cluding extremely "red" flat spectrum radio quasars (FSRQs), with 
Inverse Compton (IC) peak at keV-MeV energies, or "blue" BL 
Lacs with synchrotron peak at these frequencies, and provides a 
way to test the validity of the blazar sequence, looking for outliers 
not predicted by the sequence itself (Giommi et al. 2007). In par- 
ticular, by selecting hard X-ray luminous FSRQ at high redshift, 
i.e. with the highest intrinsic bolometric luminosities, it is possi- 
ble to collect samples of supermassive black holes (SMBHs) in the 



early Universe, thus introducing important observational constraint 
on the number density of heavy black holes at high redshift and 
hence on their formation processes and timing. 

IGR J225 17+2217 was first reported by Krivonos et al. (2007) 
as an unidentified object detected by Integral-IBIS (150 ks of total 
IBIS exposure). Swift follow-up observations were used to asso- 
ciate the source with MG3 J225 155+2217 (Bassani et al. 2007). 
MG3 J225155+2217 has been optically identified as a QSO by 
Falco et al. (1998) in a redshift survey of 177 FSRQs, on the ba- 
sis of Siv, Ly«, Cn and Cry emission lines. IGR 122517+2217 
is the highest redshift (z — 3.668) blazar detected in the fourth 
/nregra/-IBIS hard X-rays catalog (Bird et al. 2010). The source 
has a Swift-BAtf counterpart (SWIFT J2252.0+2218) in the 3 years 
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BAT survey (Ajello et al. 2009) and is present in the multifrequency 
"Roma-BZCAT" catalog of blazars (Massaro et al. 2009). 

Using XRT, IBIS and archival optical/IR data, Bassani et al. 
(2007) constructed a non simultaneous SED of the source, show- 
ing an extremely bright X-ray emission with respect to the optical 
emission (atox < 0.75), and suggested that IGR J22517+2217 
could be a rare FSRQ with synchrotron peak at X-ray frequencies, 
or a more canonical FSRQ, i.e. with the synchrotron peak at radio- 
mm frequencies and IC peak at MeV-GeV energies, but with an 
exceptionally strong Compton dominance. 

This "controversial" blazar has been studied also by Maraschi 
et al. (2008). They reanalyzed the existing Swift (XRT and BAT) 
and Integral-IBIS data, and propose a "standard leptonic one-zone 
emission model" (Ghisellini & Tavecchio 2009, see Sect. 3) with 
the peak of the synchrotron component at microwave/radio fre- 
quencies, and a high luminosity external Compton (EC) compo- 
nent peaking in hard X-rays to reproduce the SED of the source. 
This model ruled out both a synchrotron and a synchrotron self- 
Compton (SSC) interpretation for the X-ray emission. 

Ghisellini et al. (2010) included IGR J22517+2217 in their 
sample of 10 X-ray selected blazar at z > 2: the intent of the 
paper was to characterize the physical properties of these powerful 
sources, and to confirm the capability of the hard X-ray selection 
in finding extreme blazars with massive SMBH, powerful jets and 
luminous accretion disks. IGR J22517+2217 is the highest redshift 
FSRQ in their sample and shows the highest total jet power (Pj e t = 
1.5 x 10 48 erg s _1 ). 

All these previous studies have been performed through the 
analysis of the "average" X-ray spectra obtained with the INTE- 
GRAL and Swift-BKT surveys, without taking into account any 
possible flux variation of the source during the period of monitor- 
ing (5 years for INTEGRAL and BAT). In this paper we present 
the discovery of strong flaring activity, in X-ray (IBIS and BAT) 
archival data, of this extremely bright and peculiar FSRQ, and the 
modelization of both its flaring and quiescent SEDs. New Suzaku 
and Fermi data are used for characterizing the quiescent state. Our 
goal is to investigate the evolution of the SED and obtain infor- 
mation on the physical condition of the source in the two different 
states. 

The paper is organized as follows. In §2 we report the mul- 
tiwavelength data analysis of the instruments involved in the SED 
building. In §3 we describe the model adopted to reproduce the 
broad band SED, while in §4 we discuss the SED fitting of both 
the flaring and quiescent state. The summary of our results is pre- 
sented in §5. Throughout the paper, a A — CDM cosmology with 
H = 71 km s" 1 Mpc -1 , fi A = 0.73, and fi m = 0.27 is adopted. 



2 DATA REDUCTION 
2.1 New Suzaku data 

Suzaku observed the source on 2009 Nov 26 (ID 704060010, PI A. 
De Rosa), for a net exposure of ~40 ks, with both the X-ray Imag- 
ing Spectrometer (XIS; Koyama et al. 2007) and the Hard X-ray 
Detector (HXD; Takahashi et al. 2007). The XIS instrument has 3 
operating detectors at the time of the observation: the front illumi- 
nated XIS and XIS 3, sensitive in the 0.5-10 keV band, and the 
back-illuminated XIS 1, extending the low energy range to 0.2 keV. 
The HXD is instead composed of GSO scintillators and silicon PIN 
diodes. The PIN detectors observe in the 12-60 keV energy band, 
while the GSO ones can observe up to 600 keV. Data reduction 



and processing were performed using HEASOFT ?;6.9 and Suzaku 
pools vl6. Only XIS and PIN data have been used in this analysis, 
since the source is below the sensitivity limit of the GSO scintilla- 
tors. The cleaned event files produced by the data processing with 
the standard selection criteria were used. 

XIS source events were extracted from a region of radius 200 
arcsec, and the background events from an annulus (external radius 
400 arcsec) outside the source region. The response matrix and ef- 
fective area were calculated for each detector using Suzaku ftools 
tasks xisrmfgen and xissimarfgen. Given that XIS and XIS 3 have 
similar responses, their event files were summed. Only data in the 
0.5-8 keV band were considered, and the spectral counts were re- 
binned using at least 30 counts per bin, in order to allow the use of 
X 2 statistic. 

PIN data were extracted from the HXD cleaned event files af- 
ter standard screening. The tuned background model supplied by 
the Suzaku team (Fukazawa et al. 2009), was used for the "Non X- 
ray Background" (NXB) events. The background light curve was 
corrected for the lOx oversampling rate. The source spectra were 
corrected for deadtime using hxddtcor. We estimate the cosmic X- 
ray background (CXB) contribution to the PIN background using 
the model given in Gruber et al. (1999), which is folded with the 
PIN response to estimate the CXB rate. 

The XIS and PIN data were fitted with a simple power-law, 
modified by neutral absorption at the source redshift, plus galactic 
absorption, fixed to the value measured by Kalberla et al. (2005) at 
the source coordinate^ (Nh.Go.1. — 5 x 10 20 cm -2 ). The best fit 
values are reported in Table 1. The XIS and PIN data show a flat 
spectrum, with F = 1.5 ±0.1 andT = 1.5 ±0.8 respectively, with 
the PIN power law normalization being ~ 1.2 times the XIS one, a 
known cross calibration issue0. The XIS data show some curvature 
below 1 keV, that can be reproduced either by an intrinsic column 
density of Nh = (1.5 ± 1.1) x 10 22 cm -2 at the source redshift 
(compatible with the value found in Bassani et al. using XRT data) 
or with a broken power-law with break energy of 0.8 ± 0.1 keV 
and AT = 0.8. Both models give a comparable reduced \ 2 of 1.03 
and 1.04, respectively, thus the quality of the data do not allow 
do disentangle the two possibilities. In Sect. 4 we will attempt to 
disentangle these two different models using broadband data. 

The hard X-ray flux obtained by fitting the PIN data is a factor 
of ~10 (4) lower that the one reported in literature from Integral- 
IBIS (Swift-BAF) survey data (Bassani et al. 2007; Maraschi et al. 
2008; Ghisellini et al. 2010). This suggests that the source was 
in a much less active state during the 2009 Suzaku observation, 
compared to previous measurements. This lead us to reanalyze the 
archival IBIS and BAT data in order to check for the presence of 
variability in the hard X-ray source flux. In Figure 1 top panel we 
show the Suzaku XIS (red, black) and PIN (green) unfolded spec- 
tra and model best fit found for IGR J225 17+2217. For comparison 
also the IBIS spectrum (blue) is shown, to emphasize the different 
flux state. 



2.2 Integral-FBIS and Swift-BAF survey data 

IBIS data are taken from the Fourth IBIS/ISGRI Soft 7-ray Sur- 
vey Catalog (Bird et al. 2010), covering the period Feb 2003 - 
Apr 2008. Images from the ISGRI detector (Lebrun et al. 2003) for 



1 http://heasarc.nasa.gov/cgi-bin/Tools/w3nh/w3nh.pl 

2 http://heasarc.gsfs.nasa.gov/docs/suzaku/analysis/abc sec 5.4 
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Figure 1. a) Top panel: Suzaku XIS (red, black) and PIN (green) unfolded 
spectrum and model of IGR J225 17+2217. For comparison the IBIS cat4 
spectrum is shown (blue), b) Bottom panel: Hard X-ray light curve of IGR 
J22517+2217. Red squares (black circles) represent BAT (IBIS) 15-55 keV 
flux. The blue empty (green filled) triangles show the time of XRT/UVOT 
(Suzaku) observation. The cyan solid lines represent the time intervals used 
in the BAT spectra extraction. The orange dashed line represents the period 
of observation of Fermi/LAT. 



each pointing have been generated using scientific Analysis Soft- 
ware (Goldwurm et al. 2003) OS A version 7.0. Five primary energy 
bands (20-40, 30-60, 20-100, 17-30, and 18-60 keV) were used to 
maximize the detection sensitivity for sources with various energy 
spectra (for details see Bird et al. 2010). 

The BAT analysis results presented in this paper were derived 
with all the available data during the time interval 2005 Jan 19 - 
2010 Apr 14. The 15-55 keV spectra and light curve were extracted 
following the recipes presented in Ajello et al. (2008, 2009). The 
spectra are constructed by weighted averaging of the source spectra 
extracted from short exposures (e.g., 300s) and are accurate to the 
mCrab level. The reader is referred to Ajello et al. (2009) for more 
details. 

The total IBIS and BAT spectra can be reproduced with a sim- 
ple power law having T = 1.6 ± 0.6 and 1.6 ± 0.5, respectively; 



the 15-55 keV flux is Fis-sskeV = (2.5 ± 0.9) x 10" 11 and 
Fi5_55koV = (2.1 ±0.8) x 10 -11 erg cm -2 s _1 , respectively. 

The bottom panel of Fig. 1 shows the BAT (red squares) and 
IBIS (black diamonds) historical light curve of IGR J22517+2217 
from 2003 Dec 04 to 2010 Mar 23. The 15-55 keV IBIS light curve 
was extracted from the ISDC Science Products Archivqj adopt- 
ing a nominal bin size of 100 days. We converted the observed 
IBIS and BAT count rates into 15-55 keV observed flux, using 
the WebPimms HEASARC toofl assuming as underling model a 
power law with photon index T = 1.6, consistent with the values 
observed in the BAT and IBIS spectra, and assuming a constant 
cross-calibration between IBIS and BAT, equal to one. 
Fig. 1 (bottom panel) shows that: 

- the source displays quite strong long term variability in hard 
X-rays; 

- a strong flare episode occurred around Jan 2005, and the 
source reached a 15-55 keV flux maximum of (8 ± 2) x 10" 11 erg 
cm" 2 s -1 (a factor of 20 higher than the flux measured by Suzaku- 
PIN in 2009); 

- after the flare, the source faded into a quiescent state, reach- 
ing a flux that is at or below the detection limit of both BAT and 
IBIS instruments. As can be seen, the IBIS light curve is completely 
dominated by the flare, i.e. the source flux is below the IBIS detec- 
tion limit after MJD 53550 and the total spectrum extracted from 
the entire period can be considered representative of the flare. 

As a result of its different pointing strategy, BAT has a much 
more regular and extended coverage of the source, and we were 
able to characterize the source in both states. We extracted a BAT 
spectrum from the period around the flare of 2005 (2004 Dec 1 1- 
2005 Mar 21) and also one from the remaining quiescent period 
(2005 May 10-2009 Jun 18, solid cyan lines in the bottom panel of 
Fig. 1). 

The BAT flux relative to the flare state is the highest one mea- 
sured in the hard X-ray energy range (Fi 5 „s 5 k c v = (3.7 ± 0.8) x 
10" 11 erg cm" 2 s" 1 ). During this state the source is detected up 
to 200 keV, and the spectrum is characterized by photon index of 
1.5±0.5. In the quiescent state the source is detected, with signif- 
icance ~ 3<r, only up to ~ 75 keV, and the spectrum has a flux 
Fi5-55kcV a factor of ~ 15 lower than the flaring one, while the 
photon index is 1.7±1.1. Considering the large uncertainties on V 
we can consider the spectra (in flaring and quiescent state) compa- 
rable (see Table 1 for IBIS and BAT spectral analysis results). 



2.3 Fermi/LAT data 

Fermi Large Area Telescope (LAT; Atwood et al. 2009) data were 
collected from Aug 2008 (MJD 54679) to Aug 2010 (MJD 55409). 
During this period, the Fermi/LAT instrument operated mostly in 
survey mode, scanning the entire 7-ray sky every 3 hours. The 
analysis was performed with the ScienceTools software package 
version \>9rl7p0, which is available from the Fermi Science Sup- 
port Center. Only events having a high probability of being photons 
- those in the "diffuse class" - were used. 

The energy range used was 100 MeV - 100 GeV, and the max- 
imum zenith angle value was 105° . We adopted the unbinned likeli- 
hood analysis, using the standard P6_V3_DIFFUSE response func- 
tions. The diffuse isotropic background used is isotropic Jem_v02 



3 http ://w w w. isdc .unige . ch/heavens_ webapp/integral 

4 http://heasarc.nasa.gov/Tools/w3pimms.html 
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Date 


Inst. 


Exp. 


A'„ 


r 


F2-10 


logL2-10 


F15-55 


logLl5_ 55 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 



2003 Dec 04-2007 Nov 17 


IBIS 


191 




1.6±0.6 






25.1±8.9 


48.1 


2004 Dec 1 1-2005 Mar 21 


BAT Flare 


160 




1.5±0.5 






36.9±8.1 


48.3 


2005 May 10-2009 Jun 18 


BAT quiesc. 


1610 




1.7±U 






2.6±1.6 


47.3 


2007 May 01 


XRT 


40 


2.0 ± 1.5 


1.4±0.1 


2.4±0.4 


47.1 






2009 Nov 01 


XIS 


40 


1.5 ±1.1 


1.5±0.1 


1.2±0.1 


46.7 






2009 Nov 01 


PIN 


40 




1.5±0.8 






3.8±1.8 


47.4 



Table 1. Best fit model parameters for the different observations of IGR 22517+2217 analyzed in this work. The broadband continuum is reproduced with a 
simple power-law, modified by intrinsic absorption at the source redshift, plus galactic absorption Njj.Gal = 5 X 10 20 cm -2 . Column: (1) Date; (2) Instrument; 
(3) Effective exposure (ks); (4) Column density (in 10 22 cm -2 units); (5) Photon index; (6) 2-10 keV observed flux (in 10~ 12 erg s — 1 cm -2 units); (7) Log 
2-10 keV deabsorbed luminosity (erg s _1 ); (8) 15-55 keV flux (in 10~ 12 erg s~ 1 cm~ 2 units); (9) Log 15-55 keV luminosity (erg s" 1 ). 



and the galactic diffuse emission model used is glljem 
We considered a region of interest (Rol) of 15° from the IGR 
J225 17+22 17 position. All sources from the 1FGL catalog (Abdo 
et al. 2010) within the Rol of the source were included in the fit, 
with their photon indexes and the integral fluxes free to vary, plus 
a point source with power law spectrum at the IGR J225 17+22 17 
position, having photon index fixed to 2.2 (a value typical of FSRQ 
in the GeV band). The normalization of the background was also 
left free to vary. We also repeated the analysis using the source list 
of the 2FGL catalog (Ackermann et al. 2011), obtaining consistent 
results. 

IGRJ22517+2217 is located- 6° from 3C 454.3, a bright, ex- 
tremely variable source in the 7-ray sky (Ackermann et al. 2010). 
It contributes for more than 90% to the total counts in the Rol. We 
tried to exclude the period of flaring activity of 3c 454.3 from the 
data set, in order to minimize the contamination by this source. The 
results obtained for IGR J22517+2217, however, did not change 
significantly. 

IGR J225 17+22 17 is not detected in the 2 year observation, 
and the computed Test Statistic (TS, Mattox et al. 1996) is TS ~ 
3.5 in the full band. We therefore calculated the 95% upper limits in 
5 energy bands (i.e. 0.1-0.3, 0.3-1, 1-3, 3-10, 10-100 GeV), using 
the profile likelihood method. The upper-limits were corrected for 
attenuation due to extragalactic background light (EBL) through 
7 — 7 interactions (Chen et al. 2004, Razzaque et al. 2009), al- 
though only the 10-100 GeV band was found to be affected by sig- 
nificant attenuation. As a consequence of being close to 3C 454.3, 
the background around the position of IGR J225 17+2217 is higher 
than in a typical extragalactic field, making it difficult to constrain 
more strongly the upper limits at GeV energies. The upper limits 
are plotted in Figure 2, with all other data discussed in this paper. 

Given that Fermi data were collected starting from Aug 2008, 
they completely fall in the quiescent period of the source, and so 
it is not surprising that IGR J225 17+2217 was not detected. There- 
fore, they will be used to characterize the state of the source in the 
quiescent SED only, while we do not have available 7-ray data for 
the flaring state SED. 



2.4 Archival Data 

We collected archival radio and optical data from NASA/IPAC EX- 
TRAGALACTIC DATABASE |f| for the source. Radio data at 1.4 
4.8 and 8.4 GHz comes from NRAO/VLA. Optical data in J, H, 
and K bands are taken with the UFTI instrument at UKIRT (Kuhn 
2004). The UV data are taken from Ghisellini et al. (2010), where 
they corrected the Swift-UVCfT observed magnitudes for the ab- 
sorption of neutral hydrogen in intervening Lyman— a absorption 
systems. We also reanalyzed the archival XRT spectrum (average of 
4 contiguous observations), obtaining consistent results with those 
found in Bassani et al. (2007) using the same set of data, and with 
the results reported in §2. 1 for the XIS spectra. The only difference 
is in the XRT observed 2-10 keV flux, that is a factor 2 higher than 
the flux measured by XIS. 



3 THE SED MODEL 

The model adopted to fit the SED is a leptonic, one-zone syn- 
chrotron and inverse Compton model, fully discussed in Ghisellini 
& Tavecchio (2009). The assumptions can be summarized as fol- 
lows: 

- The emitting region is assumed to be spherical (with ra- 
dius R) and at a distance .Rdiss (dissipation radius) from the cen- 
tral black hole. The emitting electrons are injected at a rate Q(7) 
[cm s~ ] for a finite time equal to the light crossing time R/c. 
The adopted function Q(7) is a smoothly broken power law with a 
break at 7b and slopes si and S2 below and above 7b, respectively. 
The emitting region is moving with a velocity f3 c corresponding to 
a bulk Lorentz factor V. We observe the source at the viewing angle 

e v . 

- The external radiation sources taken into account are: the 
broad line region (BLR) photons, assumed to re-emit 10% of the 
accretion disk luminosity from a shell-like distribution of clouds 
located at a distance i?BLR = 10 17 L^/% 5 cm (Kaspi et al. 2005), 
where Ld is the disk luminosity; the IR emission from a dusty torus 
located at a distance R m = 2.5 x 10 18 Z/^4 5 cm (Elitzur 2006); 
the direct emission from the accretion disk, including its X-ray 
corona. Also the starlight contribution from the inner region of the 



http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html 6 http://ned.ipac.caltech.edu/ 
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Figure 2. a) Top panel: Spectral energy distribution of IGR J22517+2217. Gray circles and arrows represent archival radio/optical/UV data from NED. 
Empty red triangles and magenta squares represent XIS and XIS 3 data, empty green circles and orange pentagons represent PIN and BAT quiescent data 
respectively, while black arrows are Fermi upper limits in 5 bands. Filled violet squares represent XRT data, filled cyan diamonds and blue pentagons represent 
IBIS and BAT flare data respectively. The solid cyan and orange curves are the results of the modeling of the quiescent and flaring states, respectively. With gray 
lines we show the different components of the non-thermal emission: synchrotron (dotted), synchrotron self-Compton (long dashed) and external Compton 
(dot-dashed). The black dashed line corresponds to the thermal emission of the disk, the IR torus and the X-ray disk corona. The model does not account for 
radio emission, produced from much larger regions of the jet. a) Bottom panel: Zoom in the X-ray energy range for the two SEDs. Symbols as in top panel. 
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State i? diss iV nj -B r 7 b 7 max si s- 2 

(1) (2) (3) (4) (5) (6) (7) (8) (9) 



Low 570(1900) 0.045 1.06 16 70 2e3 -1 4 
High 990(3300) 0.30 0.61 16* 70* 2e3* -1* 4* 



Table 2. Input parameters of the SED fitting for the low and high state of IGR J22517+2217. Column: (1) State ; (2) dissipation radius in units of 10 cm 
and, in parenthesis, in units of Schwarzschild radii; (3) intrinsic injected power (10 45 erg s _1 ) in the form of relativistic electrons; (4) magnetic field intensity 
(Gauss); (5) bulk Lorentz factor at Rdiss', (6) and (7) break and maximum random Lorentz factors of the injected electrons; (8) and (9) slopes of the injected 
electron distribution [Q(7)] below and above 7 b ; * fixed values from low to high state SED. 



host galaxy and the cosmic background radiation are taken into ac- 
count, but these photon sources are unimportant in the case of IGR 
J225 17+2217. 

- The accretion disk is a standard "Shakura & Syunyaev" 
(1973) disk, emitting as a blackbody at each radius. The maximum 
temperature (T max ), i.e. the peak of the disk luminosity, is assumed 
to occur at ~ 5 Schwarzschild radii (Rs). Thus from the position 
of the peak of the disk luminosity, and the total luminosity of the 
accretion disk (Ld) it is possible to derive Mbh and M, once a 
value for the efficiency r\ is assumed {r\ = 0.08 for a Schwarzschild 
black hole). See Ghisellini et al. (2010a) for a discussion on caveats 
of this black hole mass estimate. 

We can estimate the black hole mass Mbh and the accretion 
luminosity L d of IGR J225 17+2217 using optical and UV data and 
upper-limits corrected for the absorption of neutral hydrogen in in- 
tervening Lyman a systems along the line of sight (see Ghisellini 
et al. 2010a for details). Given the uncertainties in the amount of 
intervening Lya systems and the paucity of data the results must 
be considered an approximation. We find Mbh = 10 9 Mq and 
L d = 6.8 x 10 46 erg s _1 . These values correspond to a disk radi- 
ating at 45% of the Eddington level. 

The BLR are located at Rblr = 8 x 10 17 cm and the IR 
emitting torus at i?iR = 2.5 x 10 19 cm. The total X-ray corona 
luminosity is assumed to emit 30% of L d . Its spectral shape is as- 
sumed to be oc v~ x exp( — hu/ 150 keV). 



4 DISCUSSION 

Our analysis shows that the extremely Compton dominated FSRQ 
IGR J225 17+2217, experienced a strong flare in the high energy 
hump in Jan 2005, and then faded in a quiescent state. In order to 
investigate the physical properties of the source we built two SEDs 
for the two different states, and fitted the data with the leptonic, 
one-zone synchrotron and inverse Compton model described in §3. 

In order to build the SED for the quiescent state of IGR 
J225 17+2217, we used the X-ray data from Suzaku (XIS and 
PIN), the BAT spectrum extracted from the quiescent period and 
FermifLKT 24 months upper limits. The archival non-simultaneous 
optical/UV data were added to our data. 

For the flaring SED we used hard X-ray data from Integral- 
IBIS and the Swift-BKT spectrum extracted from the Jan 2005 flare. 
We do not have soft X-ray data available for the flaring period. 
However, in order to put some constraint at soft X-ray frequencies 
we chose to include the XRT spectrum in the flaring SED. This 
has a factor 2 higher normalization with respect to the XIS data. 
We stress that, as observed in other samples of bright, red blazars 
(Ghisellini et al. 2010a), the flux variability is larger at higher ener- 



gies (i.e. in hard X-rays) but modest at few keV. This support our 
choice of include the XRT data in the flaring SED. 

All data points and SED model components of the flaring and 
quiescent SED are plotted in Fig. 2, top panel. Fig. 2, bottom panel, 
shows a zoom in the region of X-ray data. The strong and very 
bright hard X-ray spectrum, together with the upper limits in the 
Fermi/LAT energy range, constrain the peak of the high energy 
hump of the quiescent SED to be located at ~ 10 20 -10 21 Hz, and 
thus the form of energy spectrum of radiating electrons. The same 
energy spectrum is then assumed in the flaring SED, for which no 
MeV/GeV data are available. The corresponding synchrotron peak 
falls in a region of the spectrum localized around 10 11 Hz for both 
SEDs. 

However, in our single-zone leptonic model, the synchrotron 
and the high energy humps are produced by the same population of 
electrons. Furthermore, if the high energy emission is given by the 
external Compton process, the energies of the electrons emitting at 
the ~MeV peak are rather modest, implying a corresponding low 
frequency synchrotron peak. We then require that the synchrotron 
component peaks at low energies, close to the self-absorption fre- 
quency, and furthermore require that the thin synchrotron emission 
has a steep spectrum, whose flux is smaller than the optical archival 
data (characterized instead by a rather hard slope, that we interpret 
as emission from the accretion disk). 

Given the relative paucity of the observational data, the choice 
of the model parameters is not unique, and some further assumption 
has to be made. We indeed fix the viewing angle 6 V to 3°, to be 
close to the 9 V ~ 1/T condition. 

Note that we model both states of the source by changing the 
minimum number of parameters, given that fitting two data sets 
with the same model, allows to constraint better the model param- 
eters. In particular, we assume that the accretion luminosity does 
not change for the two states, and require also the bulk Lorentz fac- 
tor and the parameters of the distribution of the injected electrons 
(break and maximum random Lorentz factors and slopes below and 
above 7 b ) to be the same for the two SEDs. Thus, the parameters 
that are left free to vary from one state to the other are the dissipa- 
tion radius i?di ss , the injected power P( ni and the magnetic field B, 
that is proportional to l/i? diss (the Poynting flux is assumed to be 
constant). 

The results of our modeling are shown in Fig. 2, where we 
show the total flux together with the contributions of the non- 
thermal (synchrotron, self Compton, external Compton) and ther- 
mal (accretion disk, IR torus, X-ray corona) components. As can 
be seen, the curvature around 1 keV observed in the XRT and 
XIS spectra can be well reproduced by an EC component chang- 
ing (softening) slope from ~ 10 17 to ~ 10 19 Hz, disfavoring the 
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intrinsic obscuration scenario for the shape of the X-ray emission 
of IGR J225 17+22 17. 

Table 2 reports the parameters of the SED fitting for the quies- 
cent and flaring states. The main difference between the two SEDs 
is the power P ; ' n j injected in the source in the form of relativistic 
electrons, that changes by a factor ~7. The increase of P/ n j ac- 
counts for the enhanced X-ray flux in the high state. The other 
difference is the location of the emitting region R<u ss , becoming 
larger for the high state. This is dictated by the detailed modeling 
of the slope of the soft to hard X-ray spectrum, requiring an elec- 
tron distribution with a break at low energies, and another break at 
somewhat larger energies. This is accounted for, in our modeling, 
by requiring that electron of very low energies (corresponding to 
random Lorentz factors 7 < 5) do not cool in one light crossing 
time. This can be achieved if the location of the emitting region is 
slightly beyond the BLR, in a zone where the BLR radiation energy 
density is somewhat smaller. This is the reason leading to a larger 
i?diss in the high state. As a consequence of assuming a larger re- 
gion, the magnetic field is lower, following the assumption of a 
constant Poynting flux (oc B R ). The large Compton dominance 
constrains the value of the magnetic field, and in turn the relevance 
of the self Compton flux, found to be almost negligible. 

This is also in agreement with the results pointed out in Sikora 
et al. (2009): bright blazars with very hard (a x < 0.5) X-ray 
spectra and high luminosity ratio between high and low frequency 
spectral components challenge both the standard synchrotron self- 
Compton and hadronic models, while EC can easily account for 
these observed properties. 

In the analysis described above, the bulk Lorentz factor is as- 
sumed to be constant. However, the change of bulk Lorentz factor 
is often invoked to explain the variability of FSRQs. As a further 
check, we performed a new fit of both low and high states, leaving 
T as a free parameter, in addition to 7i dias and P/ nj , although the fit 
is limited by the poor number of data points 

Both SEDs are well reproduced with these new parameters 
for the low (high) state: T=15 (20); P d i ss =1700 (3500) R s and 
Log(P ; ' nj )= 43.48 (44.17) erg s -1 . In this case the values of P{ ni 
are slightly lower for both low and high states, and the differ- 
ence is lower (a factor of ~ 5), while the difference in Pdiss is 
slightly larger: from 1700 to 3500 R,s instead of 1900 and 3300 
R,s . Thus the change in T has the main effect of slightly decreasing 
the required variation of the total injected power to account for the 
observed variability, but no other substantial differences are intro- 
duced. This new fit also gives an idea of the degree of degeneracy 
in the fit parameters, due to the incompleteness of the data set, es- 
pecially in the 7-ray band. 

In Table 3 we report the logarithm of the jet power in the form 
of radiation, Poynting flux and bulk motion of electrons and pro- 
tons, in erg s~\ calculated for both SEDs. They have been calcu- 
lated from 

Pi = ivR 2 T 2 cU- (1) 

where U[ is the energy density of interest, calculated in the co- 
moving frame (see e.g. Celotti & Ghisellini 2008). As discussed in 
Ghisellini et al (201 1), the power P r dissipated by the jet to produce 
the radiation we see is almost model-independent, since it depends 
only on the observed luminosity and on the bulk Lorentz factor T. 
The power dissipated in other forms, on the other hand, depends 
on the amount of electrons (P e ), protons (P p ), and magnetic field 
(Pb) carried by the jet, which have been estimated by applying our 
specific model. Furthermore the power carried in the bulk motion 
of protons requires a knowledge of how many protons are there, per 



State 


logP r 


logP B 


logP 


logPp 


(1) 


(2) 


(3) 


(4) 


(5) 


Low 


46.04 


45.54 


45.00 


47.56 


High 


46.83 


45.54 


46.17 


48.41 



Table 3. Input parameters of the SED model fitting both the low and high 
state of IGR J22517+2217. Column: (1) State ; (2)-(5) logarithm of the 
jet power in the form of radiation (P r ), Poynting flux (Pb), bulk motion of 
electrons (P c ) and protons (P p , assuming one proton per emitting electron). 
Powers are in erg s~ 1 . 



emitting lepton. The values given in Table 3 assume one proton for 
one emitting lepton. 

From Table 3 we can see that the power P r changes from ~ 

0. 15 x Ld in the quiescent state, to P r ~ Ld in the flaring state, 

1. e. the jet requires a power comparable to the disk luminosity to 
produce the radiation we see, in the latter case. 

Tanaka et al. (2011) reported a similar behavior, based on 
Fermi-LAJ data, for the strong 2010 GeV flare of the blazar 
4C+21.35: assuming similar efficiencies for the accretion disk and 
the jet, they estimated a jet intrinsic power Lj et changing from 
~ O.lLacc in the quiescent state to ~ lL acc (beeing L acc the in- 
trinsic accretion power). They argued that these results, combined 
with the findings of Abdo et al. (2010b) on several FSRQ detected 
by Fermi-LKT, suggest a scenario in which the observed 7-ray vari- 
ability of blazars is due to the different power of the jet, that nor- 
mally represents only a small fraction of the accretion power, while 
during major flares is capable of carry away almost all the avail- 
able accretion power. Our findings are quantitatively similar, and 
thus result in agreement with this view. 

The P r however is a lower limit to the total jet power. The total 
jet power, dominated by the bulk motion of protons associated to 
emitting electrons, is P jct = P B + P c + P P = 3.6 x 10 47 and 
2.6 x 10 4S erg s^ 1 , in the low and high state, respectively. Pj e t is 
dominated by P p : if there is indeed one proton per emitting lepton, 
then the jet power is from 3 to 30 times more powerful than the 
accretion luminosity. 

It has been proposed that jets in luminous blazars may well be 
numerically dominated by pairs, but being still dynamically dom- 
inated by protons (Sikora & Madejski 2000, Kataoka et al. 2008). 
We computed the jet power due to protons assuming an upper limit 
for this ratio of one proton per 20 pairs. Above this limit the jet 
is too "light" and the Compton drag produces significant deceler- 
ation (Ghisellini & Tavecchio 2010). The one-to-one ratio can be 
assumed as a lower limit for this ratio. 

The lower limits for P p obtained in this way are 1.8 x 10 46 
and 1.3 x 10 47 erg s _1 for the low and high state, respectively. 
Therefore with this assumption, the jet power in electron, protons 
and magnetic field, becomes comparable with the radiation power. 
This translates into a total jet power of Pj ct = 2.2 x 10 46 and 
1.4 x 10 47 erg s _1 , respectively. 

The values obtained assuming one proton per lepton are ex- 
treme, even if compared with the distribution of Pj ot and Ld 
computed, with the same assumptions, for a sample of high red- 
shift Fermi/LAS and BAT blazars in Ghisellini et al. (2011, Fig. 
9), that are in the range P jct ~ 10 46 - 2 x 10 4S and L d ~ 
8 x 10 45 - 2 x 10 47 erg s _1 . To better clarify the remarkable 
behavior of IGR J225 17+2217, we note that similar values of jet 
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power have been achieved during the exceptional flare of 3C 454.3 
in December 2009 (Ackermann et al. 2010; Bonnoli et al. 2011). 
Therefore, IGR J225 17+22 17 represents one of the "monsters" of 
the high-z Universe. 



5 SUMMARY AND CONCLUSION 

Thanks to a new Suzaku observation in the X-ray energy band, the 
Fermi upper limits in the 0.1-100 GeV band, the flux selected spec- 
tra obtained through a re-analysis of the IBIS and BAT hard X-ray 
data, and other optical and radio archival data sets, we were able to 
identify a strong flare episode in the high redshift, hard X-ray se- 
lected blazar IGR J22517+2217, which occurred in Jan 2005, fol- 
lowed by a period of quiescence, extending up to the present days. 
To model the overall SEDs of the source in the flare and quies- 
cent states, we adopted a leptonic, one-zone synchrotron and in- 
verse Compton model. The optical/UV emission is interpreted as 
thermal emission from the accretion disk, plus IC from the corona, 
and reprocessed emission from the BLR. 

The curvature observed in the X-ray spectra was proposed to 
be due to intrinsic, moderate absorption (Nh ~ 2 x 10 22 cm -2 ). 
However, in the context of the broad band SED modelization pro- 
posed in this paper, it appears to be inherently accounted for by an 
intrinsic softening of the EC component around ~ 10 18 Hz. 

In both states a very strong Compton dominance is observed, 
with the high energy hump (produced by EC), that is at least two or- 
ders of magnitudes higher than the low energy (synchrotron) one. 
The high energy peak flux varies by a factor 10 between the two 
states, while the high energy peak frequency remain almost con- 
stant, between 10 20 — 10 22 Hz. The observed large Compton dom- 
inance constrains the value of the magnetic field, and hence the rel- 
evance of the self-Compton component, that is found to be negligi- 
ble in both states. The model can explain the observed variability as 
a variation of the total number of emitting electrons (a variation of 
factor ~ 7) and as a change in the dissipation radius, moving from 
within to outside the broad line region as the luminosity increases. 

In the flaring state, the jet power lower limit, represented by 
the radiative component P r , requires a power comparable to the 
disk luminosity to produce the observed radiation. The total jet 
power upper limit, dominated by the bulk motion of protons, and 
estimated assuming one proton per electron, is more than ~ 30 
times more powerful than the accretion luminosity (2.6 x 10 48 erg 
s -1 ). Such extreme values have been derived only recently for a 
handful of extreme, high redshift, hard-X/soft-7 ray selected FS- 
RQs, showing similar strong Compton dominance, and comparable 
with the value achieved by 3C 454.3 during the 2009 exceptional 
flare. 
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